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Abstract. We give a widely self-contained introduction to the mathematical 
theory of the Anderson model. After defining the Anderson model and deter- 
mining its almost sure spectrum, we prove localization properties of the model. 
Here we discuss spectral as well as dynamical localization and provide proofs 
based on the fractional moments (or Aizenman-Molchanov) method. 

We also discuss, in less self-contained form, the extension of the fractional 
moment method to the continuum Anderson model. Finally, we mention major 
open problems. 

These notes are based on several lecture series which the author gave at 
the Kochi School on Random Schrodinger Operators, November 26-28, 2009, 
the Arizona School of Analysis and Applications, March 15-19, 2010 and the 
Summer School on Mathematical Physics, Sogang University, July 20-23, 2010. 



1. Introduction 

In 1958 the physicist P. W. Anderson introduced the model which is now named 
after him to explain the quantum mechanical effects of disorder, as present in mate- 
rials such as alloys and amorphous media [9]. The most famous phenomena which 
arise in the context of this model are Anderson localization, i.e. the suppression of 
electron transport due to disorder, and the Anderson transition in three-dimensional 
disordered media which predicts the existence of a mobility edge separating energy 
regions of localized states from an extended states region. Anderson localization 
has important consequences throughout physics, in theory and experiment. Ander- 
son's work, and that of N. F. Mott and J. H. van Vleck, won the 1977 physics Nobel 
prize "for their fundamental theoretical investigations of the electronic structure of 
magnetic and disordered systems'^- 

Mathematically rigorous studies of the Anderson Model and other models of 
random operators started in the 1970s, with the first proof of Anderson localization 
for a related one-dimensional model provided by I. Goldsheid, S. Molchanov and 
L. Pastur in 1977 [23], followed several years later by a proof of localization for the 
actual Anderson model by H. Kunz and B. Souillard |29], also initially for dimension 
one. Since then the study of random operators has become an important field of 
mathematical physics, which has led to a tremendous amount of research activity 
and many mathematical results. 

While the Anderson transition and extended states are still an open mathemat- 
ical challenge, by now a good rigorous understanding of Anderson localization has 
been achieved. Several powerful methods have been found to prove Anderson local- 
ization. Important differences exist between one-dimensional and multi-dimensional 
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models, where different pliysical mecfianisms are responsible for localization effects. 
In these notes we will focus on methods which allow to prove Anderson localization 
in arbitrary dimension. Two such methods are available: The method of multiscale 
analaysis (MSA) developed in 1983 by Frohlich and Spencer [20 , and the fractional 
moments method (FMM) introduced by Aizenman and Molchanov in 1993 [5]. 

MSA has produced results in situations which are out of reach for an approach 
through the FMM, see Section l9T] for some related discussion. However, the FMM 
is mathematically more elementary, in particular for the case of the classical dis- 
crete Anderson model which will be our main focus here. Also, under suitable 
assumptions, the FMM allows to prove stronger results on dynamical localization 
than can be obtained by MSA. Therefore, in these lectures, after an introduction to 
the Anderson model and its basic spectral properties, we will discuss how to prove 
Anderson localization based on the FMM. 

After more than 50 years of physical research and more than 30 years of math- 
ematical work a vast literature with results on Anderson localization and, more 
generally, the physics of disordered quantum mechanical systems, is available. In 
these introductory lectures we ignore most of the literature as it can not be our goal 
to provide a comprehensive survey, not even of the mathematical research which 
has been done. Some book length presentations, or parts of such, which provide 
very good further reading and many more references are [TSl [TH |331 1331 [57] . 

It is assumed below that the reader is familiar with measure and integration 
theory as presented in j37j . with basic probabilistic concepts such as independence, 
and with the foundations of the theory of linear operators in Hilbcrt spaces, up to 
the spectral theorem for self-adjoint operators and consequences such as spectral 
types (absolutely continuous, singular continuous and pure point spectrum) and the 
abstract solution of the time-dependent Schrodinger equation via Stone's theorem, 
e.g. 45^ or f35|. Otherwise, we have tried to keep these notes mostly self-contained. 
For much of the first seven sections we provide full proofs. 

We do not aim at the most general known results, but rather want to demonstrate 
that simple and natural mathematical ideas can be used to rigorously establish 
Anderson localization. Many further developments of the ideas discussed here can 
be found in the literature. The references provided below can serve as a starting 
point for further reading. An ideal source for continued reading and learning the 
state of the art of much what is discussed here will be the upcoming book [8 by 
M. Aizenman and S. Warzel. 

In Section [5] we introduce the Anderson model and, as a warm-up, prove its first 
important property, namely that its spectrum is almost surely deterministic. The 
rest of these notes exclusively deals with the phenomenon of Anderson localization. 
Section [3] introduces the concepts of spectral localization and dynamical localization, 
followed by a discussion of what is known on the physics level of rigor. 

In Sections S] and [S] we prove localization in the large disorder regime of the 
Anderson model. This is done via the fractional moments method, by first proving 
in Section|4]tliat fractional moments of Green's function decay exponentially, and by 
then showing in Section[5]that this implies dynamical as well as spectral localization. 
In these sections we use methods which were developed in some of the first papers 
on the fractional moments method, e.g. [5] and |24] . In particular, these methods 
work directly for the Anderson model in infinite volume. 
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Subsequently, other methods were introduced in the literature, e.g. [T] or [6], 
which use finite volume restrictions of the Anderson model. A central concept 
here are so-called eigenfunction correlators. These methods have proven to be very 
powerful in further-reaching work, for example in dealing with the continuum An- 
derson model or mult i- particle Anderson models. Thus we introduce this approach 
in Section |6] and Appendix [X] and show how they yield an alternative proof of 
localization. 

Section [7] discusses the second main regime in which multi-dimensional localiza- 
tion has been established rigorously, the band edge regime. Among the new ideas 
needed here are the phenomenon of Lifshits tails of the integrated density of states 
near spectral edges and a geometric decoupling method to control correlations in 
Green's function. Parts of this section have the character of an outline, referring 
to the literature for some of the results used. 

Entirely written in form of an outline is Section [51 in which we discuss the ex- 
tension of the FMM to continuum Anderson models, as accomplished in ^2" and 
[12] . This requires considerable technical effort and we merely point out the diffi- 
culties which had to be overcome and mention some of the tools which allowed to 
accomplish this. 

The Anderson model and, more generally, the quantum mechanics of disordered 
media, provides many difficult future challenges for mathematicians. We discuss 
some of them in our concluding Section [9l 

Acknowledgement: The author's knowledge of random operators and, in par- 
ticular, of the Anderson model has benefitted from many other mathematicians, 
through their works as well as through personal contact. We apologize for not be- 
ing able to properly give credit for the origins of all of these benefits. But we need 
and want to make an exception for the contributions of Michael Aizenman, who 
was the driving force in the development of the fractional moments method and has 
influenced the author's way of thinking about random operators in multiple ways. 
Much of what we have to say here is based on ideas of Michael and his collaborators. 
In particular, special thanks are due to Michael and Simone Warzel for letting the 
author use some preliminary material from in Section [6] and Appendix lAl below. 

Thanks are also due to a referee for useful suggestions which improved our pre- 
sentation. 

Finally, the author would like to thank the organizers of the Kochi School on 
Random Schrodinger Operators, the Arizona School of Analysis and Applications 
and the Summer School on Mathematical Physics at Sogang University for their 
invitations. Without these opportunities to lecture on the material covered here 
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2. The Anderson Model 

2.1. The Discrete Laplacian. Below we will introduce the Anderson model as a 
discrete Schrodinger operator, acting, for dimension d > 1, on the Hilbert space 



with inner product {u,v) — X^n '"(")^("')- 
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The usual negative Laplacian — A = — J2j d'^/d^Xj is replaced by its discrete 
analogue ho, which acts on u G £^(Z'^) by 

(1) {hou){n) = - u{n + k), neZ'^, 

keZ'', |fc|=i 

where |fc| = j^i | + . . . + is the graph distance on Z'^. More appropriately, a 
finite difference approximation of — A would be given by ho + 2(1, but we neglect the 
mathematically trivial shift by 2d (which should still be kept in mind for physical 
interpretations). In physics, the Hamiltonian ho (or its negative) most frequently 
arises more directly, i.e. not as a discretization of a differential operator, in the 
context of the so-called tight-binding approximation. 

Just as the continuum Laplacian, the discrete Laplacian is unitarily equivalent 
to a multiplication operator via Fourier transform. Here we consider the Fourier 
transform 

F : L'^{[0,2Trf) ^ f{Z'^), 
which is the unitary operator given by 

{Fg){n) = {2^)-''/^ ( g{x)e-^^-- dx, 

with inverse 

(F-iu)(2;) = l.i.m. (27r)''*/2 ^ u(n)e"-". 

n£Z'i,\n\<N 

Here x ■ n = xiUi + . . . XdUd and l.i.m. denotes the limit iV — >■ oo in £^(Z''). 
A calculation shows that 

d 

F-^hoF = -2^cos(a;j), 

i=i 

where the right-hand side is understood as a multiplication operator on L^{[Q, 2tt)''-) 
in the variable x = (xi, . . . , Xd)- The function g{x) = —2 J2j cos{xj) is real-valued 
and bounded. Thus ho is bounded and self-adjoint (which also can be checked 
directly from the definition ^ without use of the Fourier transform). The range 
of g gives the spectrum of ho, 

(2) a{ho) = [-2d,2d]. 

With a bit more effort one can show that inverse images of Lebesgue-nuUsets in 
M under the function g are Lebesgue-nuUsets in M'*. Thus the spectrum of ho is 
purely absolutely continuous. 

Another similarity of ho with the continuum Laplacian is that is has plane waves 
as generalized eigenfunctions. To see this, let x e [0, 2tt)'^ and set 

(3) 0,(n):=e™-^ 
While (1)^ ^ ho acts on it via (P) as 

fe|=i 



^2cos(x-,) (j}x{n). 
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Thus (j)x is a bounded generalized eigenfunction of ho to the spectral value 

2.2. The Anderson Model. Let oj = {u)n)nei''' be a set of independent, identically 

distributed (i.i.d.) real-valued random variables indexed by n e I/^. Recall that 
this means the following, where we denote probabilities by P: 

• The (w„) are identically distributed, i.e. there exists a Borel probability 
measure /x on M such that, for all n gZ'^ and Borel sets ^ c M, 

P(w„ eA) = n{A). 

• The {ojn) are independent. Thus, for each finite subset {ni, . . . ,ne} of Z'^ 
and arbitrary Borel sets Ai, . . . , C M., 

e 

¥{iJm G Ai,...,LJne & Ae) = Yl ^i'^rij € Aj) 

£ 

i=i 

It is sometimes useful to think of a concrete way in which i.i.d. random variables 
can be realized as measurable functions on a probability space {Cl,A,¥). The 
standard construction is the infinite product space 

{n,A,F) = (g)(M,eM,/x), 

with A and P denoting the tr-algebra and measure generated by the pre-measure 

induced by /i on the Borel cylinder sets in f7 = . This is consistent with the 
notation ui = {iOn)n&<^ the components a;„ of w € are now i.i.d. random 
variables on O with common distribution /it. 

It is also convenient to work on a complete probability space ($1, A, P), which in 
the above realization is achieved by completing the product algebra A under P, for 
which the same notation will be kept. 

The Anderson Model is a random Hamiltonian /i^^ on ^^(Z''), defined for w e O 

by 

(4) {hoju){n) = {hou){n) + uJnu{n), n e Z'^. 

Introducing the random potential : Z** — >• M by T4j(n) = w„, we may also 
write 

= hQ + K;. 

Note here that is not a single operator, but rather an operator- valued function 
on a probability space. It's operator theoretic properties will generally depend on 
w. Our goal will typically be that a certain property of hi^ holds alm,ost surely or 
with probability one, meaning that it holds for w G fig, a measurable subset of Q, 
with P(no) = 1. It lies within the nature of random operator theory that the most 
interesting properties will only hold almost surely rather than for all u G Q. 

One may think of the Anderson model hi^ as the Hamiltonian governing the 
quantum mechanical motion of a single electron in a discretized alloy-type random 
medium. In this view the random potential VL,(n) = a;„, n G Z'', represents a 
solid formed by nuclei located at the sites n of the lattice Z'' and carrying random 
electrical charges a;„. Assuming that /i„ is self-adjoint, the dynamics of the electron 
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is given through the unitary group e~ " , defined via the spectral theorem, which 
provides the solution = e~^'^'^'-"ipQ of the time-dependent Schrodinger equation 
huii^it) = i4''{t), "0(0) = ipo- The possible energies of the electron are given by the 
spectrum a{hi^) of the Anderson Hamiltonian hi^. 

In the above discussion we have assumed self-adjointness of the Hamiltonian, 
which is the first mathematical fact to be checked. This is particularly easy for 
discrete Schrodinger operators such as because the discrete Laplacian ho is 
bounded and self-adjoint. 

Theorem 1. For every ui £ , the operator h^ is self-adjoint on 
D{V^) = {ue £2(2<i) . ^ |c^„u(n)|2 < oo}, 

n 

the domain of the maximal multiplication operator by the potential Vuj . 

Proof. Self-adjointness of the maximal multiplication operator by a real-valued 
function is a standard fact. Perturbation of the self-adjoint maximal multiplication 
operator by the bounded self-adjoint operator ho preserves self-adjointness with 
same domain, e.g. [45]. □ 

If we assume, as will be done later, that the distribution of the uin has bounded 
support, i.e. that 

supp /I {t e M : - e,t + e)) > for all e > 0} 

is bounded in M, then the potential is bounded and therefore defines a bounded 
multiplication operator. Thus h^i is a bounded self-adjoint operator on £^(2"^) as 
well. On the other hand, if supp/z is unbounded, then it is not hard to see that VL,, 
and thus /i^j, is almost surely unbounded. 

2.3. The spectrum of the Anderson model. Our next goal is to determine 
the spectrum of h^ . It follows as a consequence of the general theory of so-called 
ergodic operators (e.g. |14)). of which the Anderson model is a special case, that 
cr{hu) is almost surely deterministic, i.e. there exists a closed subset S of R such 
that 

o'ihuj) = S almost surely. 
Rather than proving this within the general theory of ergodic operators we will 
give a direct proof of the following result, which explicitly describes the almost sure 
spectrum of the Anderson model: 

Theorem 2. The spectrum of the Anderson model is almost surely given by 

(5) cr(/icj) = cr{ho) + supp^. 

Here the sum of two subsets A and S of R is defined hy A + B := {a + h : a E 
A,b £ B}. In particular, this means that the almost sure spectrum of h^ is a union 
of intervals, namely of translates of [—2d, 2d] by the points in supp/x. If supp /i 
doesn't have large gaps, then the almost sure spectrum of /i^ is a single interval. 

Proof. We begin with the easy part of the proof, namely that 

(6) (j{h^) C [-2d, 2d] + supp ^ 
almost surely. 

We first argue that a{V^) — {a;„ : n G 1/^} C supp/i almost surely. In fact, as 
/Lt(supp/i) — 1, for fixed n £ Z'^, ojn £ supp /it holds almost surely, i.e. on a set 
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ri„ C ri with P(r2„) = l. The countable hitersection il' of the fin also has measure 
one and for w e fi' we have cr(Kj) C supp as supp is closed. 

By a general fact from spectral theory, easily proven using a Neumann series 
argument, a bounded self-adjoint perturbation B does not shift the spectrum of a 
self-adjoint operator A by more than i.e. a{A + B) C <t{A) + 

Thus dS]) holds for lu eO.' hy 1^. 

The proof of 

(7) [-2d, 2d] + supp ^ C a{h^) 

with probability one is more involved and falls into a probabilistic part and a 
spectral theoretic part. 

For the probabilistic part, start with fixed t e supp fi, e > and N ^ N and let 

^t,N,e '■— {uj E V, : There exists a cube Ajv C Z'' of side length N such that 
\ujn ~ t\ < e for all n e A^v} • 
As t e supp/i, we have p :— fi{{t — e,t + e)) > 0. Thus, for each fixed cube A of 
side length in Z*^, the probability that |w„ — <| < e for all ri G A is > 0. We 
can cover Z'' by infinitely many disjoint cubes of side length N, where these events 
are independent. It follows that F{^lt,N,s) = 1- 

Next, let flt^e '■= f^Nefi^t.N.s- Thus P(f2t_e) = 1 and we will prove the following 
below: For each w G flt,e and a G [~2d, 2d] it holds that 

(8) [a + t~e,a + t + e]na{h^)y^9. 

Assuming that ([8]) is true, we proceed as follows: Define fit n^gN^^t.i/^, such 
that P(f^t) = 1. For w G f^t we have by dH) that, for aU a G [-2d, 2d], 

{a + t - + t + ^) n cr{h^) 9 

for all G N. As a{h^) is closed, this implies that a + t G a-{hi^) for all a G [—2d, 2d], 
and thus [-2d, 2d] + t C <j{hi^). 

For one last argument involving countable intersections of full measure sets, let 
i? be a countable subset of supp /i which is dense in supp /i and let D,q :— Clti^B^t- 
Then P(r2o) = 1 and for uj E Hq we have [—2d, 2d] + i? C a{h^). Using again that 
<7{h^) is closed completes the proof of ([7]). 

We still need to show which is the spectral-theoretic part of the proof. Let 
w G Vlt.e- Thus, by assumption, for each A^ G N there exists a cube Ajv of side 
length A^ such that |aj„ —t\<e for all n G Ajy. 

To a G [—2d, 2d] pick x = {xi, . . . , Xd) such that a — —2 cos(a;j) and consider 
the corresponding generalized eigenfunction (f>x from ([3]). Then ip^ := xAn4>x has 
finite support and, in particular, lies in £'^{Z,'^). We claim that 

(9) hmsup —— < e. 

To find a norm bound for {h^ - {a + t))^!;^ = {ho - a)xAK4>x + {V^j -t)xAN4'x, we 
first note that by assumption || (K; — t)xAjv'/'a: II ^ sllV^Afll- Moreover, as {ho — a)(px = 
0, it follows that ((/iq — ci)XAN4'x)in) is non-zero only for n close to the boundary 
of Ajv, where its values are bounded by a constant independent of N. Thus 

Wih^ - (a + i))V^^|| < CN^'-^y^ + s\\M\- 

On the other hand wc have llV'Jvll = N'^^^. This proves (|y|). 
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We conclude by a standard argument: If hu: — {a + t) is invertible, then, by ([9]), 

[|(/,^_(a + t))-i||>i. 

£ 

This imphes ^ by using the fact that for general self-adjoint operators A it holds 
that 



dist(z,cr(A)) ■ 

□ 



3. Localization Properties 

We will be interested in localization properties of the Anderson model, which 
can be described either by spectral properties or by dynamical properties of the 
Hamiltonian. 

To be more precise, let / C M be an open interval. We say that /i^ exhibits 
spectral localization in / if h^^ almost surely has pure point spectrum in J, i.e. 
/ does not contain any continuous spectrum of h^j, and its eigcnfunctions to all 
eigenvalues in / decay exponentially. 

If / is a non-trivial interval contained in the almost sure spectrum of h^^, which 
is a union of intervals, then spectral localization in / necessarily means that the 
spectrum consists of a dense set of eigenvalues (whose closure fills all of /). This 
phenomenon is very different and much more subtle than the appearance of dis- 
crete isolated eigenvalues, which is the classical situation encountered in atomic 
or molecular hamiltonians. In fact, the possibility of dense pure point spectrum 
historically can be considered as the biggest mathematical surprise provided by the 
investigation of the Anderson model. 

On the other hand, we say that exhibits dynamical localization in / if there 
exist constants C < oo and /i > such that 

(10) ¥.(snY>\{e,,e-''''-xi{K)ek)\ \ <Ce-^l^-^-|, 

forallj,fc € l'^. Here{ej}jgzd is the canonical orthonormal basis in Z'^, ej{k) = 5jk, 
and E(-) denotes the expectation with respect to the probability measure P, meaning 
E(X) = Jj^ X cflP for random variables X on $7. Both, e"'*'''^ as well as xii^uj), are 
defined via the functional calculus for self-adjoint operators. By xi we denote the 
characteristic function of J, so that Xii^i^) the spectral projection for h^^ onto /. 

Dynamical localization in the form (ITO)) is a strong form of asserting that so- 
lutions of the time-dependent Schrodinger equation huiip{t) = idtip{t) are staying 
localized in space, uniformly for all times, and thus shows the absence of quantum 
transport. Let us illustrate this by showing that dynamical localization implies that 
all moments of the position operator are bounded in time, i.e. for all p > and all 
finitely supported ip £ £'^{'Z'^), 

(11) sup|||A:|''e~'*''-x/(/ic.)V'll < oo almost surely. 
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where the position operator \X\ is defined by (|X|0)(n) = \n\(f){n). To see how (fTTj) 
follows from ([T0|) . assume that ip{k) — for |fc| > R. Then 



J2 {ej,e-'"'-xiihMm 

\k\<R 
j \k\<R 

where the last step used the Cauchy-Schwarz inequality. We can drop the square 
from \{ej, e~'*''"x/(/ii^)efe)p (as this number is bounded by 1) and then take expec- 
tations to get 



< E E |jf^E(sup|(e„e-*''-x/(Mefe)|)ll^f 

< ^^E E IjI'^^-'^'-'^m' 

3 \k\<R 

< OO. 

This implies the almost sure statement in (1111) (with square at the norm and there- 
fore also without). 

Dynamical localization is not only the physically more interesting statement 
than spectral localization (as physicists usually have little patience and limited 
appreciation for spectral theory), it is also the mathematically stronger property: 
We will show later that dynamical localization in / implies spectral localization in 
/. 

Let us discuss situations in which localization, spectral or dynamical, is expected 
physically. For this it will help to introduce an additional disorder parameter A > 
in the Anderson model and define 

(12) h^,x = ho + XV^, 

with V^{n) — LOn as above. Formally, this fits into the same framework as ([4]), using 
the re-scaled distribution 

(13) F{XuJr,eB)=^ix{B)■■=f^iB/X) 

of the i.i.d. random variables Acj„. The distribution is spread out over larger 
supports for larger A, corresponding to a wider range of possible random charges 
in an alloy- type medium. Thus A >> 1 is the case of large disorder and A << 1 
represents small disorder. 

Physicists know all of the following: 

In dimension d = I the entire spectrum of hi^^\ is localized for any value of the 
disorder A > 0. 

In dimension d > 2 the entire spectrum is localized at large disorder, i.e. for 
A » 1. 
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For small disorder A different behavior arises in dimensions d — 2 and d — 3. For 
d — 2 one still has localization of the entire spectrum, but possibly in a weaker form 
than for d = I, e.g. a small amount (or weak type) of quantum transport might 
be possible. On the other hand, in dimension d = 3 one observes the Anderson 
transition. There are localized regions near the band edges of the almost sure 
spectrum, separated by mobility edges from a region of extended states in the center 
of the spectrum. Extended states are interpreted as the existence of quantum 
transport in the sense that the moments (jlip should be infinite for sufficiently 
large p. The physical expectation for d = 3 is that this starts at p = 1/2, which 
corresponds to the presence of diffusive motion. 

Mathematically, localization has been proven for three different regimes: (i) 
for all energies and arbitrary disorder in d = 1, (ii) in any dimension and for all 
energies at sufficiently large disorder, and (iii) near band edges of the spectrum in 
any dimension and for arbitrary disorder. 

The mechanisms which cause localization in the Anderson model are fundamen- 
tally different for the one-dimensional and multi-dimensional case, which is also 
reflected in the mathematical methods which have been used to prove this. In 
d = 1 strong tools from the theory of one-dimensional dynamical systems are avail- 
able, in particular results on the asymptotics of products of independent random 
variables which allow to prove positivity of Lyapunov exponents. Large parts of the 
books [Ml and [34l are devoted to the presentation of the one-dimensional theory. 
A complete presentation of the Kunz-Souillard proof of localization for the one- 
dimensional Anderson model can be found in |18j . For a somewhat later survey of 
results on one-dimensional localization see f44l. 

As discussed in the introduction, we will focus here on methods which allow to 
prove multi-dimensional localization and, among the two methods which have been 
shown to accomplish this, focus on the fractional moments method. Using this 
method we will give a detailed proof of large disorder localization and also explain 
how it works to show band edge localization, in each case in arbitrary dimension. 

We will not discuss localization proofs via multiscale analysis. Excellent intro- 
ductions to this method can be found in [27 and [43^, while the state of the art of 
what can be obtained from Frohlich-Spencer-type multiscale analysis is presented in 
[2T] and the review [28] . We also mention the recent powerful extension of the ideas 
behind multiscale analysis in I13j . which allow to prove localization for continuum 
Anderson models (see Section [5]) with discretely distributed random couplings, a 
result which is beyond what can be obtained by the fractional moments method. 



Consider the Anderson model (jl2[) at disorder A > and in any dimension d > 1. 

Throughout the rest of these notes we will work with a stronger assumption 
on the distribution fi of the random parameters w„, namely that /i is absolutely 
continuous with density p. 



i.e. p is bounded and has compact support. In particular, this means that the 
Anderson hamiltonian h^i^x is a bounded self-adjoint operator in £'^{Z'^). 



4. Localization at large disorder 



(14) 
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Introduce the Green function as the matrix-elements of the resolvent of hoj^x^ 

(15) G^^\{x, y; z) := {ex, [huj.x - zy^ey). 
Our first goal is to prove 

Theorem 3 ([5]). Let < s < 1. Then there exists Aq > such that for A > Aq 
there are C < oo and /i > with 

(16) E(\G^^x{x,y;z)n<Ce-^^--y^ 
uniformly in x, y e 1''- and z eC\W. 

It is the appearance of fractional moments of the form E(|-|''),0<s<l, in 
the above theorem which prompted the name "fractional moments method" for 
the circle of ideas which we want to present here. The method is also frequently 
called the "Aizenman-Molchanov method", as Aizenman and Molchanov did not 
only realize that results such as Theorem [3] hold, but that they imply spectral and 
dynamical localization. These implications will be discussed in the next section. 

The following proof of Theorem [3] closely follows the original ideas from 5j . We 
start with two lemmas, an a-priori bound on the fractional moments of Green's 
function and a decoupling lemma, which contain central ideas behind the method 
and, in increasing degree of sophistication, have been used in all subsequent devel- 
opments of the method. 

Lemma 4.1 (A priori bound). There exists a constant Ci — Gi{s,p) < oo such 
that 

(17) Ex,y{\G^,x{^,y;z)n<G,X-' 
for all x,y eU^, z e<C\ M, and A > 0. 

Here 

E^,y(...)= j J . . . p{uJx) dLUx P{i^y) dojy 

is the conditional expectation with {(^u)u£Z''\{x,y} fixed. After averaging over uj^ 
and LUy the bound in (|17p does not depend on the remaining random parameters. 
Thus we get also that 

Ei\G^,xix,y;z)n<C\X-'. 

Proof. We first prove (fT7|) for the case x = y, which demonstrates the simplicity 
of the fundamental idea underlying the FMM. For fixed x € Z"*, write ui = {uj,uJx) 
where ui is short for (i-Ou) uei.'^\{x} ■ With P^^ := {cx, ■)ex, the orthogonal projection 
onto the span of Cx, we can separate the ujx and uj dependence of h^^x as 

hu,,x = hc,,x + XuJxPe^- 

The resolvent identity yields 

(18) (/i„,A - zy^ = {hc.,x - zy' - \ojxih^,x - zy'p^yh^^x - zy\ 

Taking matrix-elements we conclude for the corresponding diagonal Green functions 
that 

(19) G^^xix,x;z) = Gc.,x{x,x;z) - XujxGc,,x{x,x; z)G^,x{x,x; z) 



or 
(20) 
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Note that the latter is weU-defined since one can easily check the Herglotz property 
IniGii,A(2:, x; z)/Iniz > of the Green function. 

The important fact is that a is a complex number which does not depend on uj^- 
Thus, writing Kx{. ■ ■) ■= J ■ ■ ■ p{(^x) d^x, we find that 



(21) 



\\P\\ 



< 



SUppp I A 



Cip,s) 
A- 



with C'{p, s) independent of A and a, and thus independent of cj, z and x. 

The proof of (fTT]) for x ^ y is based on the same idea, replacing the rank- 
one-perturbation arguments above with rank- two-perturbation arguments. Write 



W = {u},UJx,l^y), P 



Pp.. and 



AWyPe^- 



Using the resolvent identity similar to above one arrives at 



(22) 
where 



z)-'P 



A + X 



UJx 

oj. 



V 



A={P{hc,x-z)-^P)-\ 



both to be read as identities for 2 x 2-matrices in the range of P. This is a special case 
of the Krein formula which characterizes the resolvents of finite-rank perturbations 
of general self adjoint operators. For the matrix A one can check that Im^ = 
■^{A — A*) < if Imz > and Im^ > if Imz < 0. It is also independent of ui^ 
and u)y. 

Using that G^^^xix, y; z) is one of the matrix-elements of P{h^ x — z)^^P, we find 

-1 ' 



Ex^y{\G^,x{x,y;z)\') < 



A + X 



OJx 





< 



A-'^E, 



IpIIL 



-A 



V 

i^x 

w, 



A'' 



LOy 



dujx duj. 



J" 



where [— r, r] is an interval containing supp p. In the double integral we change 
variables to 

U ^ -(U}j: + UJy), V ^ -(U}x - UJy), 

which gives a Jacobian factor of 2. As {uJx,ujy) G [— r, r]^ implies {u,v) E [— r, r]^ 
we arrive at the bound 

E,J\G^,x{x,y;z)n 

mil 



< 



< 



A« 

4r||HI^ 
A'' 



■C(r, s) 



Gis,p) 



-V 
V 



-ul 



du dv 
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That the latter bound is uniform in x, y and z as well as in the random parameters 



follows from the fact that the matrix 



-V 
V 



has either positive or negative imaginary part and the following general result: 
For every s G (0, 1) and r > there exists C(r, s) < oo such that 



(23) 



f \\{B~uI)-^\\''du<C{r,s) 

J — r 



for all 2 X 2-matrices B such that either ImB > or ImB < 0. 

Let us reproduce an elementary proof of this fact, e.g. Lemma 4.1 in [25j . starting 
with the observation that, by Schur's Theorem, B may be assumed upper triangular. 
We also may assume without loss that Im B >Q. 

Thus 

^11 

622 

and 

/ _ 1 61 

(B-uI)-^ = 



(24) 



(25) 



B = 



bii—u 





(bll— m)(&22-ii) 



b22 —U 



The bound ([23|) follows if we can establish a corresponding fractional integral bound 
for the absolute value of each entry of (P5|) separately. For the diagonal entries this 
is obvious. 

We bound the upper right entry of ([25]) by 



(611 - u)(622 - u) 



< 



\bi2 



(26) 

The positive matrix 



|Im((5ii-u)(522-u))| 
1 



. Imbii+Im b22 
' |bi2| 



ImB = 



Imfcii j-b 



hhi2 Im622 



2'i 



has positive determinant, i.e. det Im_B = Im5iilm622 

2 



Im(biib22) 

\bi2\ 



|6i2p/4. We thus get 



Im 611 + Im 6; 



'22 



:^12 



21m6iilni522 1 
- 2' 



-112 



The latter allows to conclude the required integral bound for (pS)) . 



The other result needed for the proof of Theorem [3] is 



□ 



Lemma 4.2 (Decoupling Lemma). For a compactly supported and hounded density 
function p as above there exists a constant C2 < 00 such that 



(27) 

uniformly in rj, j3 C. 



I ^^!rzj^Piv)dv 



< C2 
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This can be understood as a consequence of the following two facts: (i) The two 
integrals on the left hand side of ([27|) are continuous functions of 77 and /3. As both 
of them neither vanish nor diverge, the same is true for the ratio of the integrals, 
(ii) As and jr/j become large, the left hand side of ((27)) has finite limits. This 
combines to give a uniform bound in j3 and rj. The details are left as an exercise, 
or can be found in 24 . 

We are now prepared to complete the proof of Theorem [3l 



Proof. Given the a-priori bound from Lemma |4. II we may assume y ^ x. Then 

(28) = (e„ej,) 

= {ex,{hoj,x- zy^{ht.o,x- z)ey) 

u: \u—y\ — l 

= - ^ Gu,,\{x,u;z) + {Xujy - z)G^^x{x,y;z). 

u: \u—y\ — l 

Note that Guj,\{x, y; z) is the upper left entry of the matrix on the left hand side 
of the Krein formula ((22)) . Explicitly inverting the right hand side of ((22)) we find 
that 

G^^x{x,y]z) ^ -, 

AbJy — p 

where a and (3 do not depend on ujy (and it will not matter that they depend on 
A). Using Lemma 10 the bound (X;„ \an\Y < J2n ^nd ([28]) we find 

(29) Ei\G^,xix,y;z)n = 1e( 



< — E 



= ^E{\\^y-z\'\G^^x{x,y-z)Y) 
< ^ ^ E{\G^^x{x,u-z)Y). 

u: \u—y\ — l 

If none of the lattice sites u are equal to x, then the argument can be iterated. For 
given x and y one can iterate |x — y| times, in each step picking up a factor 2dG2/X^ 
after a maximum is taken over the 2c? terms in the sums over next neighbors. This 
results in a bound 

E{\G^,x{x,y;z)n< — ^ sup E{\G^^xix,u; z)n. 



For the last term we use the a-priori bound Ci/X" provided by Lemma [01 We get 
the exponential decay in ((T6)) for A > Aq if we choose Ag such that 2^6*2/ Ag < 1. 

□ 

We conclude this section by remarking that the exponential decay bound found 
in Theorem [3] also holds for finite volume restrictions of the Anderson Hamiltonian. 
More precisely, let L G N and := [-L, L]'^ n Z''. By and G^\ we denote 
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the restriction of hui.x to £'^{Al) as well as its Green function. By the same proof 
as above one finds that, for A > Ao, 

(30) E{\GtMx,y;z)n<Ce-'^^''-'^, 

where the constants C < oo and /i > are now also uniform in L. 

Moreover, in the finite volume case the bound (j30p is uniform in z € C, allowing 
for real energy. The reason for this is that the operators are finite-dimensional 
and that any given real number E is almost surely not one of their eigenvalues, which 
implicitly follows from the above proof. In the finite volume case this also holds 
for the a- priori bound in Lemma |4. II This explains why such bounds play a role in 
the FMM similar to the role played by Wegner estimates in localization proofs via 
MSA. They demonstrate that eigenvalues are sensitive to the disorder parameters. 
A good way to think of the main idea behind the FMM is that eigenvalues are 
singularities of the resolvent which move linearly under the random parameters. 
Thus the Green function can be made integrable by taking fractional moments. 

5. From Fractional Moment Bounds to Localization 

We will now discuss methods which show that exponential decay of fractional 
moments of Green's function as shown in Theorem [3] implies spectral as well as 
dynamical localization. For the sake of stating a general result of this form we 
will absorb the disorder parameter into the random parameters uj^ (re-scaling their 
distribution as in ()13|) ). Thus we consider the Anderson Hamiltonian in its original 
form Q with single-site distribution /i satisfying (fT4|) . 

From now on we will generally leave the dependence of various quantities on the 
random variable uj implicit and write h = h^j, G = G^j, etc. 

Our goal is to prove 

Theorem 4. Let I C M be an open bounded interval. If there exist s E (0, 1), 
C < oo and fi > such that 

(31) E{\G{x, y;E + ie)\') < Ce-^l^-^l 

uniformly in E £ I and e > 0, then dynamical localization in the form JJQj) holds 
on the interval I . 

As a first consequence, by Theorem [3] this implies that at sufficiently large dis- 
order A the Anderson model is dynamically localized in the entire spectrum. In 
Section [7] below, we will also use the criterion provided by Theorem |4] to prove 
band edge localization. 

The most direct way to conclude spectral localization, i.e. pure point spectrum 
with exponentially decaying eigenfunctions, from bounds such as PT|) is by the 
Simon- Wolff method. It was developed in [42] to serve a similar purpose in the 
context of multiscale analysis, where it showed that the Green function bounds 
established in [2^ indeed imply spectral localization. A short argument, showing 
that the Simon- Wolff criterion also can be combined with ([3T1) to show spectral 
localization, is provided in [S]. 

Here we will instead discuss the proof of Theorem SJ i.e. focus on how ([31]) 
implies dynamical localization. We have two reasons for doing so: First, dynamical 
localization is the physically more relevant property. Second, as we will show at 
the end of this section, dynamical localization implies spectral localization with a 
straightforward argument using the RAGE theorem. 
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There are two substantially different arguments which prove Theorem |4l In this 
section we will present a modification of an argument provided by Graf in [24j . 
This version of the argument has recently also been used in [25 to prove dynamical 
localization for the so-called unitary Anderson model. 

The second method, via the use of eigenfunction correlators, will be discussed in 
the next section. 

Graf's argument starts with the realization that fractional moments of Green's 
functions of the Anderson model can be used to bound the second moment of 
Green's function as long as a small factor (the imaginary part of the energy) is 
introduced to control the singularities of Green's function at real energy. 

Proposition 5.1. For every s e (0, 1) there exists a constant Ci < oo only de- 
pending on s and p such that 

(32) |Imz|E,(|G(x,2;;z)n < CiE,(|G(a;, y; z)^) 

for a/Z z G C \ R and x, y G Z''. 

Here Ea; denotes averaging over ujx as in the proof of Lemma 14.11 Integrating 
over the remaining variables, we see that ([32|) also holds with E^: replaced by E. 
Our proof follows the proof of Lemma 3 in [53] almost line by line. 

Proof. As in the proof of Lemma |4. II write ui — (Cj^lOx)- Keep uj fixed and consider 
the Hamiltonian 

obtained by "wiggling the potential at x" . Its Green function will be denoted by 
G^"). Similar to ([l8l to (123 we find 



and 

'''' - i^g.t) 

1 Gui(x,y;z) 



a + Guiix,x; z) ^ Gui{x,x;z) 
For the special case x = y and a = —ReGuj{x,x; z)~-^ we get from (133]) that 



1 



lmG{x,x; z) ^ 



G'''^\x,x;z) 



< 



1 



llm zl ' 



i.e. \bnG{x, x; z) ^ \ > |Imz|. Inserting this into ([55]) gives 

f34) IImz||G(")fa vz)\^ < . \Gu{x,y;z)\ 

(34) |lmz||G [x,y,z)\ _ |^ ^ ^^(^^ ^. \GUx,x;z)\ 



2 
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On the other hand, we can bound the same expression by 
(35)|Imz||G(")(x,y;z)p < |Imz| ^ |G(")(x,2/';^)P 

- |Imz|(e,,(/i(")-z)-i(/i(")-z)-ie,) 

= |Imz|(e,, - - (M") - z)-i]e,) 

= IniG(")(x,x;z) 
_ |ImG^(a;,a;;z)~^[ 

where the last step used p3p with x — y. 

For t > one has min(l,i^) < . Using this to interpolate between ([34]) and 
([35| we get 



M ii..^iig-'(^.=/;^)p^ ir:g:i:::;:i::i. -t|gf- 

We will now use the following "re-sampling trick" , which has the effect of creating 
an additional random variable (here a) to average over. For a non-negative Borel 
function / on R, 



(37) / I f{uJx + a)p{'^x + a) da p{uJx) di^x 

f{u}^ + a)p{(jJx + a)p{(jJx) du)^ da 
f{ujj;)p{LUj;)p{uJx - a) du)x da 
f{^x)p{^x) (^J p(wj;-a)da^ dujj; 
f{ujx)p{ujr^)dujj;, 



where the integration order was interchanged in the first and third steps and trans- 
lation invariance of Lebesgue measure was used in the second. 
Choose f{uJx) = \G(ib,u]^)ix,y; z)\'^ , then ((37)) and (p6| yield 

(38) |Imz|E,(|G^(x,2/;z)|2) 

= \lmz\E^(^J \G^'''>{x,y;z)\^p{uj^ + a)da^ 

~ V ' ' Guj{x,x;z)\'' J \a + G^{x,x;zy^\'^ J' 

We now use Lemma 15.21 below with w — G^{x, x; z)^^ to conclude 

\lmz\E,{\GUx,x;z)\^)<CE,{\GUx,y;z)n 
with a constant G < oo which only depends on supp p, but not on a;, y and z. □ 

In the above proof we have used 
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Lemma 5.2. There exists a constant C = C{p) < oo such that 



uniformly in w G <C and lOx G supp p. 

Proof. Using \wY < lal** + \a + w]" , we need two estimates: 
(i) 

|Imw| / ^^'^•^ ^ J^^ da < 7r|||a|''p(a;^ + a)||oo 
J \a + w\'^ 

< 7r(Kn!p|U + |||A|V(A)|U). 

(ii) 

|Imz«| /f^^^i^da < min(— ^^,q|p|U|Imu;n 
J |a + w| |lmw|^ 

< c\\p\\l-^ 



I OO 



□ 



We now complete the proof of Theorem HI 



Proof. Consider the mixed spectral measures px,y of h, the complex Borel measures 
defined by 

(39) Px,y{B) ^ {ex,XB{h)ey) 

for Borel sets B C R. The total variation \px,y\ of Hx^y is a regular bounded Borel 
measure which can be characterized by 

(40) \l^x,y\iB) = sup / g{X)dnx,yW 

g -.R^C Borel 
lffl<l 

= sup \{ex,g{h)xB{h)ey)\, 

|9|<1 

e.g. [33. The particular choice gt{x) = e^**^ in (|40|) shows that 

(41) \^lx,y\{I)>s^lp\{ex,e-''\I{h)ey)\. 

teK 

Therefore Theorem 3] will follow from a corresponding exponential decay bound for 

e(Im..,|(/)). 

As / is an open bounded interval, it follows from Lusin's Theorem ([37 ) that one 
can replace Borel functions in (I40p by continuous functions with compact support 
in /, 

(42) ,^1(7) = sup \{ex,9ih)ey)\. 

For g £ Cc{I) it follows by elementary analysis (using that g is bounded and 
uniformly continuous) that, uniformly in A G M, 



g(A)= hm - f ^^l^dE. 

E^o+ TT J (A - Ey + £2 
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By the spectral theorem this imphes 

{ex,g{h)ey)^ ^ j 9{E){c^'-<:,{h ~ E - iey^{h - E + ie)''^ey) dE. 
This allows to estimate the expected value of (|42|) by 

< EMiminf- f K'^^' ~ ^ ' ^^^^^Miez, {h - E + ie)-^ey)\ dEj 

< liminf- /'^(E(£|(e,,(/i-i?-ie)-ie,)n)'/' 



■{E{e\{e,,{h^E + te)-'ey)\')f^ dE, 

where, in this order, Fatou, Fubini and Cauchy-Schwarz (on E) have been used. 
Now Proposition 15.11 can be applied, allowing to bound further by 



< liminf i fy"{E{\G{x,z;E + ie)f)y/'^{E{\G{z,y;E~ie)\'y^^dE 

e-s-0+ TT J J ^ — ' 

< y^^-Mlx-z|/2^-MU-d/2 

z 

In the last step the assumption of Theorem |4] was used (which also applies to 
\G{z^ y;E — ie)\ = \G{y, z;E + ie)\). The elementary bound, based on the triangle 
inequality, 

^-fj.\x-z\/2^~p,\z-y\/2 ^ ^-p,\x-z\/4:^~fj.\x-y\/i^-'p,\z-y\/4 

and another use of Cauchy-Schwarz (on the z-summation) complete the proof of 
Theorem |4l 

□ 

It deserves mentioning here that we have actually proven a stronger result than 
dynamical localization in the form (|10p. The above proof shows that for an open 
interval / on which (|3T|) holds there are constants C < oo and /i > such that 

(43) E{\fi,J{I))=E{ sup \{ex,9{h)xi{h)ey)\) <Ce-^^^-y^ 

g : H ^ C Borol 

Isl < 1 

for all x,y G Z'^. An interesting special case is g = 1, where (|43)) establishes 
exponential decay of correlations in the spectral projection xi{h)- Consequences 
of this for the conductivity of an electron gas in response to an electric field have 
been discussed in [4] . Another consequence is mentioned at the end of this section. 

Next we show that dynamical localization implies pure point spectrum via the 
RAGE- Theorem. The underlying idea is very simple: The RAGE-Theorem char- 
acterizes states in the continuous spectral subspace as scattering states (in time- 
mean). Dynamical localization excludes scattering states and thus continuous spec- 
trum. 



Proposition 5.3. Suppose that dynamical localization in the form U0\) holds in 
an open interval I . Then hi^ almost surely has pure point spectrum in I . 
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Proof. For a discrete Schrodinger operator h — ho + V in £'^{Z'^) let Pcont(h) be the 
projection onto its continuous spectral subspace. Then the RAGE-Theorem, e.g. 
[15] , says that for every -0 e 

(44) ||P,„„,(/i)x/(/i)V|p- lim hm T ^ ||;t{|,|>^^e-'*''x/(/^)V'f • 

i?,— ^OO T— i-OO Jq 1 LI I— J 

If ^ has finite support, say supp^ C {\x\ < r}, then 

IIX{|x|>fl}e-**''X/(/i)^f < IIX{|x|>fl}e"**''X/(/i)X{|.|<r}llllV'f 

< \{^..e-^''\i{h)ey)m\\ 

\x\>R,\y\<r 

where dropping a square is allowed as ||x{|x|>fl,}e~'*''x/(^)X{|2;|<r} II < 1- 
Taking expectations in (j44p implies, after using Fatou and Fubini, 

(45) mPcont{K)xi{hM'') 

^ „ f% E E(l(e.,e-*^-x/(Me,)|)||V'f. 

" |2;|>fl,|i/|<r 

By (Uni) we have E(|(ex, e-**'''-x/(/i^)ey)|) < Ce-'^l^""?^! uniformly in which 
bounds the right hand side of (|45|) by 

< lim C V e-'^l^-^l = 0. 

\x\>R,\y\<r 

We conclude that Pcont{hu>)xi{hu))''l' — for almost every uj and every -0 of finite 
support. The latter are dense in £'^{1/^) and thus Pcont(huj)xi {^ui) — almost 
surely, meaning that the spectrum in / is pure point. 

□ 

We note that the above proof of pure point spectrum docs not imply exponential 
decay of corresponding eigenfunctions. It is shown in [5^ how this follows directly 
from exponential decay of fractional moments pil) , using the Simon- Wolff- method 
It can also be deduced from by considering g{h) = fe(ft-), E £ I, using 
the result from that almost surely all eigenvalues of h^, in / are non-degenerate. 
For details on this see Section 2.5 of [5], where a corresponding argument for the 
continuum Anderson model is provided which also applies to the discrete Anderson 
model considered here. 

6. Finite Volume Methods 

The localization proof provided in Section [5] proceeds directly in infinite volume, 
i.e. does not require to consider restrictions of the hamiltonian h to finite subsets of 
Z'^. However, it is also possible to work in finite volume, prove the relevant bounds 
with volume-independent constants, and then deduce localization by taking the 
infinite volume limit (sometimes referred to as "thermodynamical limit"). This 
has conceptual advantages such as having to only deal with discrete spectra, thus 
allowing to express functions of the hamiltonian by eigenfunction expansions and, as 
described at the end of Section |4j to directly study Green's function at real energy. 
Moreover, finite volume methods have proven very robust under generalizations, 
for example in the extension to continuum Anderson models which we will discuss 
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in Section m For these reasons we will use this section to provide a different proof 
of dynamical localization, i.e. Theorem |4] above, using finite volume methods. 

Many of the ideas involved here can already be found in the Kunz-Souillard ap- 
proach to localization .29. for the one-dimensional Anderson model. They were first 
combined with the fractional moment method in 1 to prove dynamical localiza- 
tion for the multi-dimensional Anderson model. A central object are so-called finite 
volume eigenfunction correlators, arising from eigenfunction expansions. Eigenfunc- 
tion correlators are also used in similar form in proofs of dynamical localization via 
multiscale analysis, see [33] or [^5] and references therein. 

The methods to be described here are not completely disjoint from the methods of 
Section[5l As before, we consider the mixed spectral measures fj,x,y of h introduced 
in ([39]) as well as their total variation lux.yl given by (j40|) . As will become clear in 
(|49|) below, \iJ,x,y \ can be considered as an infinite volume eigenfunction correlator 
for h. We will find bounds for it by finding bounds for finite volume eigenfunction 
correlators which hold uniformly in the volume. 

Let /i^^ be the restriction of h^j to A^ = [—L, LY n and denote its Green 
function by G^^ . 

Proposition 6.1. Let < s < 1 and I an open bounded interval. Then there exists 
C ~ C {s , p, d) < oo such that 



Results of this form were first used in implicit form in [T and later stated more 
explicitly in [6^ . The exact statement given here as well as its proof below and in 
Appendix |A] follow notes provided to us in private communication by M. Aizcnman 
and S. Warzel. They used similar results also in [7 . 

Based on (|4T|) . we see that Proposition 16.11 may be applied to provide a proof of 
dynamical localization in / in situations where it can be shown that 



with constants which are uniform in L and E Cz I. In fact, dynamical localization 
follows under the somewhat weaker assumption that the energy average over / of 
the fractional moments of Green's function is exponentially decaying. However, in 
all our applications we have uniform pointwise bounds available. For example, as 
discussed at the end of Section 01 a bound of the form (|T7)) holds on the entire 
spectrum for sufficiently large disorder, thus providing a second proof of dynamical 
localization in this regime. 

Proof of Proposition \6.1l We start by reducing the claim (|46p to properties of finite- 
volume spectral measures. We again use the characterization (|^^ of \fJ.x,y{I)\ for 
open bounded intervals /. Strong resolvent convergence of to h implies for 
continuous g of compact support that {cx, g{h^'^)ey) — >■ {cx, g{h)ey) and thus, by 



(46) 




(47) 



E(|G•^(x,^/;i^;)^<Ce-^l-~^l 



(48) 




Here the regularity assumption on g can be dropped since h^'^ has discrete spec- 
trum. 
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Let be the restriction of h^^ to the reducing subspace T-Lx for h^^ generated 
by Cx and let Px be the orthogonal projection onto T-Lx- Then Cx is a cyclic vector 
for and all eigenvalues E of /i^ are simple. Thus we may label the corresponding 
normalized eigenvectors by tpE- We use the notation -0^ also for ®0 in ^^(Ai) = 

T~Lx ® "^x^ ■ 

By expanding into eigenvectors we get 

\{ex.g{h^'')ey)\ = |(e,, g(/.^)P,e,)| 
=■■ QL{x,y;I), 

and, in particular, 

sup \{ex,g{h^'')ey) \ < QL{x,y;I). 
Iffl<i 

The latter will be referred to as eigenfunction correlators. Using Fatou's lemma we 
conclude from (|48|) that 

(49) E(|^,,J(/)) <liminfE(Qi(a;,y;/)). 

In order to establish a relation to the fractional moments of Green's function we 
will also introduce fractional eigenfunction correlators through 

(50) QUx,y;I,r):^ ^ l'Al(a;)P"'l^l(y)r 

Eeina(h^) 

for < r < 2, noting that Ql{x, y; I) — Ql{x, y; 1, 1). We claim that for < s < 1, 

(51) EQL{x,y-I)<iEQLix,y;I,s))^. 

To see this, interpolate s < 1 < 2 via 1 = - + 2 with the conjugate exponents 
p = 2 — s and q — f-rf ■ Applying Holder to expectation as well as to summation 
yields 

EQLix,y;I,l) < {EQL{x,y; I , s))^ {EQLix,y; 1,2))^ . 

This implies ([51]) after observing that QL{x,y;I,2) = Y,Eeincr(h^) \i'E{y)? ^ 1- 

We will now be able to relate the fractional eigenfunction correlators to fractional 
moments of Green's function by showing that there exists a constant C — C{s, p, d) 
such that 

(52) EQL{x,y;I,s) <C jE{\G^-{x,y;E)ndE. 

This, combined with (gS]) and implies 

In the proof of ([5^ we will use the fractional eigenfunction correlators 
QL,v{x,y; I, s) which are defined as in ([50)) . but with the summation being over 
the eigenvalues and eigenfunctions of + vPe^ ■ Note that, cyclic vector 

for h^^ + '^Pe^ is the same as the restriction of h^'^ + wPe^ to T-Lx and that Cx is 
a cyclic vector for this operator for all values of v G K. This makes Proposition lA.2l 
in Appendix |X] applicable to our situation, which we will now use to finish the proof 
of Proposition [nH] by invoking a resampling argument. 
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For this note that f ^^"-^u du is continuous and non-vanishina; as a function of 
a G K. Thus there exists a constant C = C(s, p) < oo such that 

p{a) 



(53) '-^ < C for all a G M. 

f ^^"^ du ~ 

Writing uj — {(li,ujx) and denoting the expectation over a) by E, we get 

(54) E(Q^(x,y;/,s)) = E [ Qf'^-\x,y;I,s)p{iu,)du;x 

Jm 

< CE [ ( [ Qf'-^\x,y;I,s)^^^)piu)du. 
Jr\Jr \u-0Jx\'J 

After the change of variable uJx v := cj^; — u we see that the inner inte- 
gral is equal to J Q''^'^\x,y; I,s) By Proposition IA.2I this coincides with 
// l^(i^ti)(^' y'' ^)\^ Inserting into (|54)) we arrive at 

E(0^(x,2/;/,s)) < CE f j \Gfi;^^^{x,v-E)\^ dE 
= CE j^\Gt{x,y;E)\''dE. 

7. LiFSHiTS Tails and Band Edge Localization 



□ 



7.1. Band edge localization. So far the only regime in which we have proven 
localization is the case of large disorder A >> 1 in (|12p . In this section we consider 
the Anderson model in its original form (|4]), i.e. at fixed disorder. Our goal is a 
proof of localization at energies near the bottom of the spectrum. The arguments 
involved can be easily modified to show the same near the upper spectral edge. 

For notational convenience we will assume that the density p of the distribution 
of the LOx satisfies supp p = [0, Umax]- We also write Eq = —2d, which according to 
Theorem [2] becomes the bottom of the almost sure spectrum 

of h^. 

Our localization proof will again proceed via showing exponential decay of the 
fractional moments of Green's function: 

Theorem 5. For every s G (0, 1) there exist S > 0, fi > and C < oo such that 
(55) E{\GM: y;E + ie)\') < Ce"^!"-''! 

for all x,y eZ'^, E e [Eq, Eq + 6) and e > 0. 

As discussed at the end of Section SI our methods again yield a bound on the 
finite volume Green function as in (150]) , uniform in the volume and in i? G [Eq , Eq + 
S], allowing for e = 0. Thus we can conclude spectral and dynamical localization 
at the bottom of the spectrum from either of the methods discussed in Sections [S] 
of [6l working in infinite volume or in finite volume. 

Compared to the case of large disorder, we face essentially two new difficulties, 
which are illustrated by the shortcomings of Lemmas 14.11 and 14.21 The a-priori 
bound from Lemma |4T] is still valid and will be used. But, as the disorder A is fixed. 
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we can not hope that the a-priori bound also provides a "smallness mechanism" , 
which can be used to iteratively prove exponential decay. We will again proceed by 
iteration, but a different mechanism will be needed to get it started. Also, it will 
turn out that we need a different decoupling argument. Lemma |4.2[ which was used 
in this context in the proof of Theorem [3l is too case-specific and will not work for 
the geometric situations which we will encounter here. 

7.2. Lifshits tails. Physically, the new smallness mechanism is provided by the 
fact that the bottom of the spectrum Eq = inf S is a fluctuation boundary. This 
means that finite volume restrictions of h^^ rarely have eigenvalues close to 0. To 
make this precise, as before let h^^'^ be the restriction of to £'^{Al), Al = 
[-L,L]'^nZ'^. 

Lemma 7.1. For every 13 G (0, 1) there are rj > and C < oo such that 

(56) P(inf aiht^ )<Eo + L-^)< CL'^e-"^""^' 
for all LeN. 

To illustrate why this means that small eigenvalues are rare, let us assume that 
the \Al\ — (2L + l)"* eigenvalues of /if; are uniformly distributed in E. Then the 
smallest eigenvalue should be no larger than CjU^. But by ([55)1 this is extremely 
rare for large L. In fact, the methods used to prove ([56)) can also be used to prove 
that the integrated density of states N{E) of h^j satisfies Lifshits-tail asymptotics 
near the bottom of the spectrum: 

(57) 7V(i;) <Ce-''l^~^°l"'"', 

which is much "thinner" near Eq than the corresponding IDS No{E) = C\E—Eq\'^/'^ 
of the Laplacian Hq. 

For detailed proofs of Lemma 17.11 as well as (l57l) we refer to [27] or [43] , with 
the latter working in the setting of the continuum Anderson model (but applicable 
to the discrete model as well). Here we only briefly outline the reasons behind 
Lemma |7. II By the variational principle 

(58) infa(/^^ -i^o) = inf ((/l^ - £;o)^, 

= inf ({{h^'' - Eq)^,^) + u:.Mi)?)- 

Note that both terms on the right hand side of ([58| . the kinetic and potential 
energy, are non-negative. In order to find a low lying eigenvalue, they both need 
to be small. By reasons of the uncertainty principle, small kinetic energy requires 
that if) is approximately constant, ip{i) ^ C — |Al|^^/^, to be normalized. For such 

the potential energy is approximately X^igAl ^^/I^^I' which by the central limit 
theorem with large probability is close to the expected value E(cjo) > 0. The event 
^jWi/|Ai| < < E(a;o) is a large deviation and has probability exponentially 
small in jA^j. 

The weakest part of the above heuristics is the reference to the uncertainty 
principle. Slightly different ways to make this rigorous are provided in [27j and 
[43] . both requiring arguments which control the separation of the second lowest 
eigenvalue from the lowest eigenvalue. The proof provided in [27] (going back to 
work from the 1980s) uses Temple's inequality in this context, while [43] uses an 
argument based on analytic perturbation theory. 
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In the context of proving Theorem[5l Lemma mi provides a first step, a so-called 
initial length estimate: 

Lemma 7.2 (Initial Length Estimate). For every s G (0, 1) there exist C < oo and 
rj > such that 

(59) E{\G^^{x,y;E + ie)\') < CL'^e""^"'*"^" 

for all L eN, X, y € Al with \x - y\ > L/2, E e [So, Eo + iL-^/Crf+z)] and e > 0. 

Proof. Let /3 G (0,1) and, motivated by Lemma [7. 11 define the "good" and "bad" 

sets as no := {uj : inf cr(/i^^ - Eq) < 1'^} and := rj^. Then 

(60) 

n\G^' [x, y;E + ie)n - E{\G^- (x, y; E + ie)\^xna) + HIG''^ i^, E + *e)|^Xf2 J- 

Pick p > 1 sufficiently small such that sp < 1 and let q be conjugate to p, 
i + i = 1. Holder applied to the second term on the right hand side of (|60)) gives 

(61) Ei\G''-ix,y;E + ze)\'xnn) < (E(|G^M^, 2/; + ^e)!'^)) '^>(^^b)'/« 



< CL'f''e-^^'''\ 



where we have used the a-priori bound from Lemma 14.11 as well as the probability 
bound from Lemma 17.11 The first term on the right hand side of ()60p concerns 
the event where E has distance at least \L~^ from the bottom of the spectrum, 
which allows to conclude exponential decay of jG^^ (x, y; E + ie)\ in \x — y\ from a 
Combes- Thomas estimate (see e.g. [21]), giving 

(62) E{\G^^x,y;E + ie)\'xaa) < GL-^^e-^"!"-^!/^' 

< GL^^e-^i^"' 

for constants 77 > 0, ryi > and G < 00. 

The choice /3 = 2/(2 -f d) leads to equal exponents in ((6T|) and (|62l) . which 
combine to give ((59)) . □ 

7.3. Geometric decoupling. We will eventually fix L ~ Lq, choosing Lq such 
that the right hand side of (|59l) is sufficiently small (how small still to be deter- 
mined). After making this choice we will pick 6 = ^Lg thus determining 
the interval [Eq^Eq + S] in which Theorem [5] establishes localization. In order to 
derive the exponential decay bound (|55]) from this we have to develop a decoupling 
method which will allow to proceed iteratively, splitting the path from x to y into 
segments of length Lq. 

The description of this so-called geometric decoupling method will fill the remain- 
der of this section. Our argument will closely follow a construction introduced in 

El- 

In addition to h{^^ , consider h^^ , the restriction of h^: to £^(A|^), where = 
Z'^\Al. Let 

This means that 

(63) /i,, = /i(f')+T(^), 

where T^^) is the operator containing the "hopping terms" introduced by the dis- 
crete Laplacian between sites of Al and A^. We write Gi^\z) = (/lif^ — z)~^. 
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More precisely, the matrix-elements of T^^^ are 
,64) Ti^'M = {o!- It"''^'" 

where is the boundary of defined as the set of pairs (u, u') with |u — = 1 
and It G Al, m' G A^, or vice versa. 

We now perform a double decoupling of the resolvent G^j = Gc^{z) by using the 
resolvent equation for (|63|) twice, first at L and then at L + 1: 

(65) G„ = Gi"^^ -Gi'^^T^'^^G^ 

= G(f ) - G(f )T(^)G(f + Gi^)r(^)G^T(^+i)G(f +1). 

Equations of this form are often referred to as geometric resolvent identities. 

By translation invariance it suffices to prove (|55p for a: = 0. U \y\ > L + 2, 
then the first two terms on the right hand side of do not contribute to the 
matrix-element Gi^{0,y;z) and thus 

G.(0,y;z) = (eo, G^f ^T^^^G.Tf^+DGif+i^e,) 

= E E Gi^'(0,?^;^)G.(i.',i;;z)G(f'+i)(z;',y;z). 
(«,M')eri, (■u/u')ert+i 

For s £ (0, 1) we get 

(66) E(|G^(0,2;;z)n 

< ^{\G!^HO,u;z)GUu',v;z)Gt^^^'{v',y;z)\'y 
{u,u')eTL 

Here we have replaced G^^ by G^^ as and u are both in A^. Similarly, Gif^^"* 

was replaced by G^^^^ . For fixed {u, u') and {v, v') consider the corresponding term 
on the right hand side of (|66p and note that the first and last of the three factors 
are independent of and ujy. Thus, in taking the expectation we may integrate 
over uju' and ujy first and use Lemma |4. II to conclude 

(67) E (|G^ (0, u; z)G^{u', v; z)G;i^+^ {v' , y; z)|^) 

(68) < GE{\GtHQ,u;z)nGt^'+'{v\y;zW) 

(69) = GE(|G^(0,u;z)r)E(|G:^'-+^(z;',2/;z)r). 

In the last step we have used that the remaining two factors in (|68p are stochastically 
independent. Now let z = E + ie with E e [Eo,Eo + ii^^/^^+d)] rj,^^^ 

we may 

estimate the first factor in by the bound obtained in Lemma 17.21 and, after 
inserting into (|66p. find 

(70) E{\G^{0,y;z)n<CL''-'e-^^"''"^'' ^ E(|Gf-+^ (t/, y; z)^. 

\\v'\\^=L+2 

We want to use ([70]) as the first step in an iteration. The second step 

would consist in finding a bound for K{\Gi^'"^^ {v' ,y; z)\^ similar to the bound for 
E(|Gaj(0, y; z)]") given by ((70|) . with v' serving as the new origin. A problem arises 
from the fact that the underlying domain is not any longer Z'^, but A'i_^i. Iterating 
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this would result in more and more complex geometries and we would be faced with 
the problem to check if all the constants involved in the estimates leading to ((70|) 
can be chosen uniform in those geometries. 

An elegant way around this is the following result of [6], see Lemma 2.3 there, 
which allows to bound the depleted Green function G^'^^'^ in terms of the full 
Green function Guj'- 

Lemma 7.3. There exists a constant C = C{s,p) < oo such that 

E{\Gi'^+'\v',y;z)n<Ei\GUv',y;z)n + G E(|G.K, y; z)n. 

||«'||o==i+2 

The proof of this starts from the geometric resolvent identity G^f^^'' — G^ + 
^(L+i)2i(L+i)^^^ In the resulting Green function expansion over {u,u') G r^+i 
crucial use is made of the bound 

Ei\Gi'^+'Hv',u;zW\GUu',y;z)n < CE(|G.K, z)n. 

The proof of this uses another special case of the Krcin formula similar to (but 
this time tracking the dependence on all four variables ujyi, uJu, ^u' and ujy) and 
a decoupling argument similar to Lemma 14.21 (but for functions of two variables 
which are linear separately in each variable) . For the details we refer to [B] . 
Lemma 17.31 combines with ([70]) to yield 

(71) E{\G^{0,y;z)n<CL^''e~'^'^"''"^"' sup E{\G^{m,y; z)^ 

11"! ||oo<i + 2 

for some constant G < oo. With that constant we fix L = Lq such that p 
GL^''e-''^o^*''''" < L We also choose 6 := ^L^^^'-^^'^^ now. For E G [Eq^Eq + S] 
we can use (|7ip to start an iteration, 

E{\G^{ui,y;z)n < p sup E(|G„(u2, y; ^)|'), 

||«2||oo<2(-Lo+2) 

and so forth. This iteration can be carried out approximately \y\/LQ times before 
the chains ui, U2, ■ ■ .may reach y. After this number of steps we use the a-priori 
bound from Lemma l4.1l to bound the last fractional moment in the chain. We have 
proven Theorem [S] with exponential decay rate /i — \ log p\/Lq. 

8. The Continuum Anderson Model 

It took somewhat more than a decade to find a generalization of the fractional 
moment method to continuum Anderson models. Our goal in this section is to 
explain why this took so long and how it was eventually done. Here our presentation 
will be less self-contained than in previous sections. We will outline the new ideas 
which were needed and refer to the literature for details. 

The main difficulty is that the rather elementary arguments from rank-one and 
rank-two perturbation theory, which worked so well for the discrete Anderson 
model, fall far short of applying in the continuum. In the latter, each single site 
potential is a perturbation of infinite rank, which at best has certain compactness 
properties relative to the Laplacian. To make the central ideas behind the frac- 
tional moment method work in this setting required a much deeper understanding 
of some of the operator-theoretic aspects involved. Here we will follow the works [2] 
and [12], where these questions were settled. Earlier work in [55] extended certain 
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aspects of the fractional moment method to continuum models, but still relied on 
finite-rank perturbation arguments by, for example, considering continuum models 
with random point interactions. 

For our presentation here we choose to work with the deterministic background 
operator 

(72) Ho = -A + Vo 

in I/^(]R'^), where Vo is a real-valued, Z'^-periodic potential in L°°(R'^). Let Eq := 
inf cr(i?o) denote its spectral minimum. 

A continuum Anderson-type model is then given by 

(73) H^^Ho- J2 

where uj = (wn)n6Z'' is an array of i.i.d. random variables with bounded density p 

such that SUppp = [OjWmaa;]- 

The single-site potentials Un (x) — U{x~n) are translates of a non- negative bump 
function U characterized by the existence of < ri < r2 < cxd and < ci < C2 < oo 
such that 

(74) CiX{|2;|<ri} < < C2X{\x\<r2}- 

The spectrum of H^^ is almost surely deterministic, 

a{H^) = S a.s., 

and 

El := inf E = inf a(Ho - Umax ^ Un) 

n 

is characterized by choosing all couplings maximal and thus, due to our sign- 
convention, the potential minimal. It can be shown under the assumption (j74p 
that the spectral minimum is strictly decreased by the random potential: Ei < Eq. 

We will use the notation x„ = XAi(n): where Ai(ri) refers to the unit cube in R'' 
centered at n G Z'^. 

The following theorem is a special case of a result in [12]. Similar results were 
first obtained in [2,, where a "covering condition" of the form 

(75) U > cxo, c > 0, 
was required for the single-site potential. 

Theorem 6. Let d < 3 and < s < i. Then there exist 5 > 0, ^ > and C < oo 
such that 

(76) mXk{H^ -E- 2e)-'xi\n < Ce-^l'^"^! 
for all E e [El, El +6], e> and kj e . 

In Theorem [6] we use the norm of the localized resolvent Xk{H^ — E — it)~^xi 
(sometimes called a "smeared Green function") as a continuum analogue of the 
discrete Green function G{x, y;E + ie). This has also been found to be the correct 
object to consider in continuum extensions of multiscale analysis. 

Without going into the details here (which for the continuum case can be done 
similar to what was described in Section |6l see [2]), we state that exponential decay 
of fractional moments of the smeared Green function, as established in ([76]) , implies 
spectral and dynamical localization: 
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Corollary 8.1. Under the assumptions of Theorem\^ the following holds: 

(a) For almost every lo, H^j has pure point spectrum in [Ei,Ei + S] with expo- 
nentially decaying eigenfunctions. 

(b) There are constants /j, > and C < oo such that 



for all k,£ & , with the supremum taken over Borel functions 17 : M — > C. 

The overall approach to proving Theorem |6] is similar to the proof of Theorem [5] 
in the previous section. The main steps are: 

(i) A priori-bound: It can be shown that to every E2 G {Ei,Eq) and < s < 1 
there exists C < 00 such that 



uniformly in e [Si, £^2], e > and kJf^U^. 

Note here that, as opposed to the discrete case Lemma BTTl the a-priori bound 
is only shown for energies below the spectrum of the unperturbed operator Hq. 
This is a consequence of not requiring the covering condition (j75p for the single-site 
potential. If a covering condition holds, then it was shown in [2] that the bound 
([75)) holds at all energies, with a constant C on the right which grows polynomially 
in E. 

(ii) Lifshits tails: The bottom Ei of the almost sure spectrum is again a 
fluctuation boundary and close analogues to Lemmas 17. II and 17.21 as well as Lifshits 
tail asymptotics ([57]) of the IDS hold in the continuum, see e.g. [43 and [2 . As in 
the discrete case, this provides the start of an iterative procedure for the proof of 
exponential decay in (l76l) . 

(iii) Geometric decoupling: The geometric decoupling procedure described 
at the end of Section [7] can be carried out similarly in the continuum. Additional 
technical difficulties arise mostly due to the fact that the required geometric re- 
solvent identities (compare (1651) ) are less straightforward in the continuum. One 
consequence of this is the restriction of Theorem [5] to s < 1/3, which is due to the 
need of an additional three-factor Holder bound used in the decoupling procedure. 
Also, elementary decoupling bounds such as Lemma 14.21 have to be replaced by a 
more systematic construction involving resampling of the random variables aj„ near 
the surfaces at which the decoupling is carried out. For details in the setting of 
Theorem ini see [T^ . 

The only one of the above three points which we want to address in some more 
detail is the a-priori bound ([TS]). as the existence of such a bound can be seen as 
the crucial test for the possibility of using the fractional moment method in the 
continuum. 

For simplicity, we only consider the "diagonal" case k — £ — here and will 
assume the covering condition ([75]). We will discuss reasons why we could hope 
that 



(77) 




(78) 



E{\\xk{H^ - E - ie)-\ir) < C 



(79) 



supE(|l[/(ff^ ~E~ iey^UW) < 00 



e>0 
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for energies near inf E. Under the covering condition, this imphes the same result 
with U replaced by xo- 

When trying to implement ideas similar to the ones used in the proof of 
Lemma 14.11 we are faced with having to find an analogue to the Krein formula. 
It turns out that this is done by the identities known from Birman-Schwinger the- 
ory. Write 

u = (w, Wo), H^^ Hcj- ujqU. 
Then, at least formally, it is easy to derive by the resolvent identity that 
(80) U'/^{H^ - zY^U^I'' = {Abs - iool)-' 

in i^(supp[/), with the Birman-Schwinger operator 



(81) Abs - {U'^'{H^ 



It can be justified that the inverses in (jSTj) and (|80l) exist and that Abs is 
maximally dissipative. Here an operator A is called maximally dissipative if it 
is dissipative, i.e. Im(0, A^) > for all (j) in its domain, and it has no proper 
dissipative extension. This can also be characterized by the fact that {e^^^}t>o is 
a contraction semigroup. 

The identity ([80l) looks promising since the right hand side separates the de- 
pendence on clIq from the dependence on u. Indeed, if the bound ()23|) could be 
generalized from dissipative 2 x 2-matrices to general maximally dissipative opera- 
tors B, then it would immediately give us ([75]). While extends to dissipative 
N X A^-matrices, the bound C(r, s) on the right will become iV-dependent and di- 
verge for A'^ — >■ oo, as is seen by choosing B to be a diagonal matrix with entries 
1, . . . , iV. Thus it is not possible to directly extend ([23]) to the Hilbert space setting. 

However, the extension to the Hilbert space setting becomes possible if additional 
Hilbert-Schmidt multipliers are introduced. This is most naturally stated in terms 
of a closely related weak-L^-bound: 

Theorem 7. Let T-Lq and Tii be separable Hilbert spaces, let A be maximally dissi- 
pative in TLq, and let M : T-Lq Hi be a Hilbert-Schmidt operator. Then 

(a) the boundary value 

M(A -v + iOy^M* lim M(A - w + ie)-^M* 

e-!-0 

exists in Hilbert-Schmidt norm for almost every w G M, 

(b ) there exists a constant C < oo ( independent of A and M ) such that 



C\\M\\ 



1 

HS 



(82) \{v eR: \\M{A-v + iO)-Hr\\Hs >t}\< ^ 

for all t > 0. 

In ([82|) I ■ I denotes Lebesgue measure and || • \\hs the Hilbert-Schmidt norm. 

Part (a) is well known in mathematical physics and has been frequently used in 
scattering theory. Less well known is part (b), which describes the value-distribution 
of the boundary values guaranteed to exist by (a). This was proven in [33] (see also 
an appendix in [2] for a reproduction of the proof), based on the weak-L-'^-property 
of the Hilbert transform of Hilbert space- valued functions, the latter being a quite 
classical result in harmonic analysis. 
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The weak-L^-bound §^ can be turned into the fractional moment bound 



(83) / \\M{A~v + iO)-Hr\\'HsP{v)dv <C{s,p)\\M\^^' 



HS' 

where the constant C(s, p) can be chosen uniform for aU Hilbert-Schmidt operators 
M and maximaUy dissipative A. This is done by the standard layer-cake integration 
argument: If F{v) := \\M{A -v + iO)-'^M*\\Hs, then 

\Fivrpiv)dv < IIpIU / \Fiv)\^dv 

Jsuppp 



IpIIoo / suppp: \F{v)\' >t}\dt. 

JQ 



By ([5^ the integrand is bounded by mm{C,C\\M\\jjg/t^/'}, where C = | supp p|. 
Splitting the integral at the i-value where C = C\\M\\jjg/t^^^ leads to 
When trying to use (|83p for a proof of ([79]) we see from ([50]) that 

(84) - zy'U = U'/'iAss - ujolY^U^'^. 

This leaves us with one more problem to deal with: The multiplication operator 
[/^/^ in L?{W^) is not Hilbert-Schmidt. In fact, multiplication operators with non- 
vanishing functions in the continuum are never compact. 

The key to solving this last problem is that U^/"^ is relatively Hilbert-Schmidt 
with respect to —A (meaning that [/-'^/^(— A-|- 1)^^ is Hilbert-Schmidt), at least for 
d < 3, see e.g. [JS. Arguments as typical in relative perturbation theory allow to 
split the left hand side of (IM|) into a sum of terms, some of which trivially satisfy 
a fractional moment bound, while others include additional multipliers which lead 
to the Hilbert-Schmidt property required in (j83p . These arguments only work at 
energies below the spectrum of the unperturbed operator i/o, which is the reason 
for the corresponding assumption which we made when stating (|78p . 

For further details on these relative perturbation arguments as well as on the 
"off-diagonal" case =^ £ in ([75)) we refer to [T^j and conclude our sketch of the 
proof of Theorem [S] here. 

9. Open Problems and directions for future work 

To conclude this introduction into the theory of Anderson localization, we men- 
tion some open problems and discuss some wide open issues which mathematicians 
need to understand better in the future. Here we will not restrict ourselves to 
further developments of the fractional moment method, but will address broader 
aspects of the quantum mechanical description of disordered media. We will be 
relatively brief here and note that a more complete and more detailed recent dis- 
cussion of open problems in this field can be found in [10] . In particular, we do not 
attempt here to give complete references to related works. 

9.1. Singular distributions. Consider the discrete and continuous Anderson 
models and H^^, but allow for singular distributions of the random coupling 
parameters w;, i e 'L'^. The most extreme case would be the case of independent 
Bernoulli variables, i.e. P(a;i — a) — V{LLii = b) ~ 1— p. This models the physically 
interesting case of a two-component alloy. Both, the fractional moments method 
and the Frohlich-Spencer multiscale analysis, fail to provide localization proofs in 
this situation. The reason for this is that both methods to a large extend use local 



32 



GUNTER STOLZ 



averaging arguments in the random parameters, as demonstrated very clearly by 
the proof of Lemma 14.11 above. While it is possible to deal with Holder-continuous 
distributions, the Bernoulli case it out of reach for the traditional approaches. 
However, Bourgain and Kenig [13| have shown 

Theorem 8. Consider the continuum Anderson model defined by |7^[ ) and 
|73| ) with Vo = and independent Bernoulli random variables (i-Oi). Then is 
spectrally localized near Ei ~ inf S. 

Their proof is based on a substantial enhancement of the multiscale analysis 
approach and, in particular, a deeper understanding of the underlying averaging 
mechanisms (such as the role of the so-called Wegner estimates). It has also been 
shown in how the argument provided in [13' can be used to handle (a;^) with 
arbitrary non-trivial distribution. 

However, somewhat surprisingly, the same question remains open for the discrete 
Anderson model Q with Bernoulli distributed random couplings. The technical 
reason for this is that [13) uses subtle unique continuation properties of the eigen- 
functions of Schrodinger operators which are not available for lattice models. 

More generally, one can easily imagine various other models of random operators 
where the random parameters naturally have discrete distribution and where the 
available mathematical methods fail to prove localization. One such model would 
be discrete Laplacians on random subgraphs of the edges of Z'^. An open question 
is to decide if in the supercritical percolation regime, where the graph has a unique 
infinite component, the Laplacian has localized spectrum. For a recent survey on 
these models see [32] 

9.2. Extended states. Every list of open problems in random operator theory 
needs to mention the somewhat embarrassing fact that mathematicians are still 
far from understanding the physically conjectured extended states regime in the 
three-dimensional Anderson model. 

A proof of the existence of continuous (or absolutely continuous) spectrum or of 
diffusive solutions to the time-dependent Schrodinger equation for this model would 
be an important break-through. Here we would like to mention another way to 
characterize the Anderson transition from localized to extended states, namely the 
level statistics conjecture. In fact, this is how physicists can numerically distinguish 
the two regimes, which provides the most important evidence for the correctness of 
the physical heuristics explaining the transition. 

According to the level statistics conjecture it is possible to distinguish the lo- 
calized and delocalized regimes by considering the statistical distribution of the 
eigenvalues (viewed as point processes) of finite volume restrictions of the Ander- 
son model. Localized states should be characterized by Poisson statistics of the 
eigenvalues, while in spectral regions with extended states the finite volume eigen- 
values should show GOE statistics. The latter it a special kind of level repulsion 
observed for Gaussian orthogonal ensembles in random matrix theory. 

In the spectral regimes where mathematicians can establish localization, it has 
also been verified that the finite volume eigenvalues are Poisson distributed, see 
[31], [30], [17] and [22]. However, regarding GOE statistics in the Anderson model, 
as little is known as for other possible characterizations of extended states. 

As discussed in the lectures by L. Erdos at this School [19], GOE statistics is 
a rather universal phenomenon observed in large classes of random matrices, e.g. 
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so-called Wigner random matrices. The most apparent difTerence between Wigner 
matrices and the Anderson model is that for the latter randomness is restricted 
to the diagonal matrix-elements while in Wigner matrices all entries are random. 
Understanding the transition between Anderson models and random matrices, for 
example by considering random band matrices with an increasing amount of off- 
diagonal random entries, could provide important insights into the localization- 
delocalization transition in the Anderson model. 

9.3. Electron-electron interactions and many-body systems. The Anderson 
models discussed above are one-electron models, which ignore interactions between 
electrons (as well as interactions between nuclei, which are considered as affixed to 
the lattice sites). Quite recently, Anderson- type models for a fixed number TV of 
interacting electrons in a random background have been shown to have localization 
properties. Chulaevsky and Suhov [151 E] have done this by an extension of mul- 
tiscale analysis, while [7] povides similar results based on the fractional moments 
approach. 

Let us give one example of a result which can be obtained by both approaches, 
where we do not try to state the most general result. An A^-particle discrete 
Anderson-type model can be defined as 

JV 

y:\y-x\ = l 3 = 1 

where e ^^(Z^'*), x = {xi, . . . ,xn) e Z^'' and y € Z^"^. As above, the {lO:,:)^^icl 
are i.i.d. random variables with bounded, compactly supported density. Assume, 
for simplicity, that U{x) is a two-particle interaction term of finite range, 

U{x)^ 5Z '^{xj—Xk)j supp<f> finite. 

l<j<*;<Ar 

Then localization holds for large disorder: 

Theorem 9. //A is sufficiently large, then hi^^ is spectrally and dynamically 
localized at all energies. 

Note that a version of dynamical localization suitable for A''-particle systems has 
to be used here, see [7]. 

While such results provide important mathematical progress, condensed matter 
physicists will object to the above model because it keeps the number of electrons 
fixed. The physically correct system to look at would be a model with positive 
electron density. For example, one could consider the restriction of hj^^ to £^(A;^), 
where A;^ = A^ x . . . x Al , and study its properties if N goes to infinity together 
with the volume, i.e. iV ^ L'' — )■ oo. This is mathematically wide open. It is not 
even clear how the concepts of localization and extended states should be defined 
in this setting, with spectral theoretic terms most likely not being the correct 
language any more. Instead dynamical descriptions will have to be used. One 
particularly interesting question, which is still challenging even to physicists, is if 
electron interaction enhances or reduces localization effects. It has been argued in 
the physics literature that such effects could be crucial to understand localization 
for two-dimensional disordered systems, the critical case. 
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Appendix A. Basic Rank-One Perturbation Theory 

Here we first provide some facts from rank-one perturbation tlieory. They can 
be considered as finite-dimensional, elementary versions of some of the much more 
profound insights behind Simon- Wolff theory, e.g. [4? or Chapter 11 to 13 of |40j . 
Then we establish a general relation between fractional moments of Green's function 
and eigenfunction correlators, which was used in Section |6] above. 

Let H be a finite-dimensional Hilbert space, dim'H = N, Hq a, self-adjoint oper- 
ator in "H, and ip a normalized cyclic vector for hg. By P — {ip, ■)ip we denote the 
orthogonal projection onto spanji^s}. Consider the family of rank-one perturbations 

hv := ho + vP, V e M, 

of ho- Then is a cyclic vector for all hy. All hv have simple eigenvalues, which 
we label as 

Ei{v) < E2{v) < ... < ENiv). 
Also consider the self-adjoint operator 

hoo := PhoP on D{hoo) = M^. 
Application of the Gram-Schmidt procedure to (p, A(p, A^ip, . . . shows that Aip — 
{ip, Aip)ip is a cyclic vector for hoc- In particular, /iqo has simple eigenvalues which 
we label as 

£;i(oo) < £2(00) < ... < En-i{oo). 
In the following we will also use the notation EqIoo) := —00, En{oo) :— 00. 

Lemma A.l. (i) For every k G {1, . . . , N}, Ef^{v) is analytic and strictly increas- 
ing as a function of v G M. 

(ii) For every w G M; the eigenvalues of hy and /iqo are intertwined as 

(85) Ei{v) < Ei{oo) < E2{v) < i;2(cx)) < . . . < EN^iiv) < EN-iioo) < En{v). 
(Hi) For every k G {1, . . . , N}, 

lim Ek{v) = Ek-i{oo), lim Ek{v) = Ek{oo). 

Proof, (i) Let ipk{v) denote normalized eigenvectors of h^ to Ek{v), k = 1, . . . , N . 
That p is a cyclic vector for means that {ipk{v), p>) 7^ for all fc G {1, . . . , N}. 
By analytic perturbation theory, the functions Ek{v) are analytic with 

(86) E'kiv) = {Mv),PM^)) = \{My),v)\^ ^0- 

(ii) We will prove this by the variational characterization of eigenvalues of hy 
and hoo, e.g. Theorem XIII. 2 in jM!, which says that for k G {1, . . . , N}, 

(87) Ek{v)^ sup inf {f,hyf) 

V cn / e v-^ 



and for fc G {1, . . . , - 1}, 



dim V = - 1 11/11 = 1 



Ek{oo) ^ sup mf {g,hoog). 

dimV = k - 1 llsll = 1 

Note that in the infimum in (|87|) the orthogonal complement V-^ is taken with 
respect to H, while V-^ in ([55)1 is taken with respect to {y'}'''. By definition of h^o 
and hy we also have that 

(89) {f,hoof)^{f,hyf) for all / G {^}^. 
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We first show that 

Ekiv)<Ekioo) for all A: e {!,..., A^}. 

This is trivial for k — N. For k < N—1, let y be a subspace of ?^ with dim V ~ fc— 1. 
Then n {ip}^ is a subspace of {</?}^ of dimension at least N — k. Thus it is 
the orthogonal complement of a subspace W of {tp}^ of dimension at most fc = 1. 
Therefore, by (p]) and (p j) . 



SfcM > inf if, hoof) = inf , ifjij) 

11/11 = 1 11/11 = 1 

> inf {f,Kf). 
f e 

11/11 = 1 

As this holds for every subspace V oiH with dim^ = fc — 1, ([87]) implies Ek{oo) > 
Ek{v). 

Next we will show that 

^fc(oo) <Ek+i{v) foraUfce{0,...,iV-l}, 

which is trivial for fc 0. Let fc > 1 and V C {</'}^ with dimV = fc* — 1. Then 
V := span{v3} ®V with dimV = fc and = {0} © V^. Thus, by dHZ]) and 

(EH), 

Ek+\{v)> inf {f,hyf)^ inf {g,hocg)- 

f ev^ g sv^ 

11/11 = 1 llsll = i 

As C {if}'^ with dim!/ = fc — 1 was arbitrary, ([88|) implies > Ek(oo). 

Strictness of all inequalities in ([85l) now is a consequence of (i). 

(iii) Here we use the following general fact, which can be proven using Schur 
complementation (see e.g. [TT] for a description of this method): For the self-adjoint 
2 X 2-block matrix 

A B 
^ B* D 

let E a{D), then 

( A + vI-EI B _ ( ^ 

|.|™oo B* D-EI ) -\{) {D- EI)-^ 

Applying this to the 2 x 2-block representation of /iq in span{(/3} {p>}^ shows that 
{K - EI)~^ -> © (/loo - EI)-^ as oo for every E ^ a{hao)- 

Using that for self-adjoint operators A, 

we conclude that for every E G cr(/ioo) there exists a function E{v) such that 
E{y) e a{hy) for all and lim^^oo E(v) ^ E. li E ^ Ek{oo) for fc = 1, . . . , - 
1, it follows from the results of (i) and (ii) that E{v) ~ Ek{v) for v sufficiently 
large, i.e. limi,_^oo i?fe(t;) — Ek{oo). Similarly, it follows that \imy^^oc> Ek+i(v) = 
Ek{oo). Ei{v) —oo as V ^ —oo and En{v) — )■ oo as f — > oo follows easily by 
minimizing/maximizing the quadratic form of hy. 

□ 
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To given % G s G (0, 1) and open interval / C M we consider the fractional 
eigenfunction correlators 

(90) Q.(v5,x;/,s):= E \{Mv),^)\'-'\{Mv),x)\'- 

k:E,{v)eI 

Proposition A. 2. The fractional eigenfunction correlators satisfy the identity 

(91) / RjikijilAdv^ [ \{^,iho-Er\)\^dE. 

Jr \w Ji 

Proof. As ip is cyclic for h^ and thus {4>k(v),f) 7^ for all k, we can rewrite ([UP)) 

as 



(92) q4p,x;I,s)= ^fe(^) 

k:Ekiv)eI 



where we have also used (|86p . Observe that 

{ho - Ek{v))Tpk{v) = {hy - Ek{v) - vP)-tJjk{v) = -'v{ip,'ipk{v))ip- 
Thus, using that Ek{v) (7{ho) for all w ^ 0, 

(^fc(«),x) = {{ho-Ek{v))M^)Aho-Ek{v))-'x) 
= -v{ipk{v),(p){(p,{ho - Ek{v)y^x)- 

This allows to further rewrite ((92|) as 

(93) ^^f:^'"^ = E i?U«)l(^,(/^o-i?.(^'))-^x)r 

k:Ek{v)£l 

This will allow to prove (j9ip by integration. Here we may assume that / C 
{Ek^i{co),Ek{oo)) for a fixed k G {1,...,A} (from which the general case fol- 
lows easily). In this case ([M)) says that 

Q,{^,xJ.s) _ f £;^(i;)|(^,(/io-Sfc(i'))-^x)h if£;fc(i^)G/, 
\ 0, else. 

Integration yields 

Qv{'P,X'-,I, s) 



I 

Jr 



dv = / Ei{v)\{ip,{ho-Ek{v))-\)\''dv 

v:Ek(v)eI 



\{ip,{ho-E)-\)YdE, 
which used the substitution v ^ E ~ E}^{v). □ 

We conclude by noting that it is possible to prove a result corresponding to 
Proposition IA.21 without assuming that the Hilbert space is finite-dimensional, as 
long as the spectral measure for ho corresponding to x is purely singular and bound- 
ary values {if, {ho — E — i0)^^x) of Green's function are used. A corresponding 
streamlining of Aizenman's original arguments in [1] has been provided by Simon 
in [41] for the unitary models considered there and will be presented for the self- 
adjoint setting in [8]. 
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